Fractional Flux Periodicity in Doped Carbon Nanotubes 
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An anomalous magnetic flux periodicity of the ground state is predicted in two-dimensional cylin- 
drical surface composed of square and honeycomb lattice. The ground state and persistent currents 
exhibit an approximate fractional period of the flux quantum for a specific Fermi energy. The period 
depends on the aspect ratio of the cylinder and on the lattice structure around the axis. We discuss 
possibility of this nontrivial periodicity in a heavily doped armchair carbon nanotube. 
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In the Aharonov-Bohm (AB) effect, the phase of wave- 
functions is modulated by a magnetic field, thereby man- 
ifesting the quantum nature of electrons. One of the di- 
rect consequences of the AB effect in solid state physics 
is the persistent current in a mesoscopic ring 1,2,3,4 . The 
persistent current is an equilibrium current driven by a 
magnetic field threading the ring, and it generally shows 
the flux periodicity of <5> = hc/e. On the other hand, 
there exist some systems where the period becomes a 
general fraction of &o, revealing an interference effect be- 
tween channels. At present, theoretical investigations of 
the fractional flux periodicity has been done in a two- 
dimensional (2D) system composed of a square lattice 
and in a one-dimensional (ID) system. Cheung et al£ 
found that a finite length cylinder with a specific aspect 
ratio exhibits the fractional flux periodicity in the per- 
sistent currents. The same configuration with an addi- 
tional magnetic field perpendicular to the cylindrical sur- 
face was analyzed by Choi et al&. They reported that a 
fractional flux periodicity appears in the persistent cur- 
rents. Its period is mainly determined by the additional 
perpendicular flux, but is also dependent on the number 
of lattice sites along the circumference. We have shown 
in the previous paper that torus geometry exhibits the 
fractional flux periodicity depending on the twist around 
the torus axis and the aspect ration. As for ID systems, 
Kusmartsev et ali& reported on fractional AB effect in 
a certain limit of the Hubbard model and Jagla et aliS 
found that correlations change the fundamental periodic- 
ity of the transmittance of an AB ring with two contacts. 
In ID, there is only one channel and therefore there exists 
no coherence effect between channels; hence the origin of 
the fractional flux periodicity differs from that of two- 
dimensional system. 

However, since these geometries have not yet been real- 
ized experimentally, it is important to examine if existing 
materials with cylindrical geometry can be used to detect 
such an interference effect. In this paper, we will show 
that an approximate fractional flux periodicity appears 
in a honeycomb lattice cylinder, which is realized in a 
single wall carbon nanotube (SWNT} 1 ™. The fractional 
periodicity requires a heavy doping, shifting the Fermi 
energy up to the energy of the transfer integral (2.9 eV). 



The persistent current in doped nanotubes was also the- 
oretically examined by Szopa et aliii, though they do not 
mention the fractional periodicity. 
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FIG. 1: (left) Geometry of a 2D cylinder in the presence of 
the AB flux, $. (right) We consider 2D cylinders composed 
of a square and an armchair nanotube. 



Here we consider the interference effect in a 2D cylinder 
(Fig. ^) composed of a honeycomb lattice, and calculate 
the ground state energy and the persistent current. The 
persistent current (I pc (N$)) is defined by differentiating 
the ground state energy (Eq(N<s>)) with respect to the 
magnetic flux penetrating through a hollow core of the 
cylinder (See Fig. UJ as I pc (N$) = -dE (N<s,)/d$. We 
show that the fraction of flux period ($>o/Z) depends on 
the aspect ratio of a cylinder: l/Z = ji/It for a doped 
(n, n) armchair SWNT with It being the number of unit 
cells along the cylindrical axis. Persistent current can be 
observed experimentally via the induced magnetic mo- 
ment of the system, which was recently demonstrated in 
a SWNT by Minot et aii£. 

First we consider a finite length 2D cylinder composed 
of a square lattice, as studied by Cheung et al&. The 
length of the cylinder and the circumferential length are 
Na and Ma, respectively, where a is the lattice constant. 
By solving the nearest-neighbor tight-binding Hamilto- 
nian with hopping integral t, we obtain the energy eigen- 
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value as 
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where 1 < n < N and h 1 < m < — . (1) 
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In the presence of the AB flux parallel to the cylinder 
axis, E nm changes according to the gauge coupling which 
can be obtained by substituting m — > m — N$, where N$ 
is the number of flux, defined by 7V$ = $/$ (see Fig.Q. 

When we consider a half- filling system (Ep = 0), the 
ground state energy Eq(N^) is given by 

E Q (N*) = J2 E En m (N$), (2) 

n=l m=[-A„+A r 4 ,] + l 

where A n = ^ (l — j^fjj and [a;] represents the largest 
integer smaller than x. From Eq. PJl. the persistent cur- 
rent is expressed as 

«^ = E E m s H M )■ 

n=l m=[-A»+iVs]+l v 7 

(3) 

For M 3> f , we can rewrite Eq. as 

Jp C (^) _ 2 ^ sinjTrJVgj 
u i=i J 

where Iq = 2ei/(Msin -p) and we ignored the correction 
of 0(M~ 2 ) to the right hand side. The coefficient, Cj is 
given by 

N /o A \ 

^=£ S in(^)cos(27rj<). (5) 

71=1 ^ ' 

When the aspect ratio satisfies M/2(N + 1) = l/Z with 
integer Z, cos(27rZ^4„) = 1 holds for all n and Cj satis- 
fies Cj+z = Cj. Furthermore, in the limit of N ^> 1, Cz 
becomes large compared with Ci,--- , Cz-i since they 
are proportional to 0(1/ N) and Cz oc O(N). Since 
only j = Z, 2Z, ■ ■ ■ are dominant in Eq. we get an 
approximate fractional ($>q/Z) periodicity of I pc when 
M/2(N + 1) = l/Z which becomes exact in the long 
length limit N 3> f , as was found by Cheung et al£ 

Let us apply this result to doped armchair SWNT. 
The carbon nanotubes can be specified by the chiral 
vector, Ch = nai + ?na 2 , and the translational vector, 
T = iiai+t 2 a 2 , where (n,m) and (ii,i 2 ) are integers and 
ai , Si2 are symmetry translational vectors on the planar 
honeycomb lattice^ with |ai| 2 = |a.2 1 2 = 2ai • a 2 . The 
chiral vector Ch specifies the circumference of the cylin- 
der, and the unit cell of the nanotube is defined by two 
mutually perpendicular vectors Ch and T. The length of 
the cylinder is specified by a vector ItT. We decompose 



the wave vector as K = //iKi + (^2/^)^-2 where (Xi and 
fj,2 are integers, and the wave vectors Ki and K2 are de- 
fined by the condition 10 : C h ■ Ki = 2tt, C h ■ K 2 = 0, 
T • K] =0 and T • K 2 = 2tt. With these definitions, 
the energy eigenvalue of the valence electrons can be ex- 
pressed by 

£7^2 = -V-k Vl + 4 cos X cos Y + 4 cos 2 Y, (6) 

where V n (= 2.9eV) is the transfer integral for nearest- 
neighbor carbon atoms and the variables X and Y are 
defined by 

X = -£- (-ft- t 2 )( Ml -j\r*) + (n-m)g), (7) 

Y=^(-(t 1+ t2)(fii-Ni,) + (n + m)^y (8) 

where N c — mt\ — ni 2 . This implies that the motions 
around the axis (fj,i) and along the axis (/i 2 ) couple with 
each other. 

Without doping (i.e., at half- filling) , the energy disper- 
sion relation of carbon nanotube has only two distinct 
Fermi points called K and K' Fermi points. In this case 
only two channels touch the Fermi level and the ground 
state energy shows the AB effect with the periodicity 
of and corresponding situation has reported ex- 

perimentalljiiSii^. We consider a special doping which 
shifts the Fermi level from to Ep — — (the case 
of Ep = Vjr can be analyzed in the same way and the 
results are the same because of the particle-hole symme- 
try). The energy dispersion relation around the Fermi 
level is approximately given by 

Bum ~E F - 2V 7r cosF(cosX + cosF) . (9) 

We note that the energy dispersion relation is similar to 
that of the square lattice except for the factor of cosy. 
Here, we adopt the armchair SWNTs (n = m) because 
we can set t\ = 1 and t 2 = — I and therefore X is pro- 
portional to only fix and Y to /x 2 , which results in that 
cos X + cos Y becomes the energy dispersion relation of 
a square lattice. 

We first show the numerical results of the ground 
state energy Eo(N<p) for the doped armchair SWNT. In 
Fig. Ufa), we plot Eq(N&) of doped armchair (30,30) 
SWNTs with different lengths. The ground state en- 
ergy shows an approximate fractional flux periodicity de- 
pending on the aspect ratio defined by lp/n. We can 
see that the approximate fractional periodicity (3>q/Z) 
corresponds to Z = lp/n. However, there also exists 
the approximate fractional flux periodicity for a shorter 
nanotube, It < nX. In Fig. Hfb), we plot the ground 
state energy for (80, 80) armchair SWNT with l T = 48 
(Z T |T| re 11.6[nm]) and l T = 96 (l T \T\ ps 23.2 [nm]). 
The approximate fractional periodicity corresponds to 
l/Z = 1/3 and 1/6, respectively, in which Ix/n is 3/5 
and 6/5. 

Now we derive the fractional flux periodicity of per- 
sistent current in doped armchair SWNTs analytically 
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between the ground state energies for the two energy dis- 
persions. Figure 13 shows the occupied states (gray color) 
for (a) the hypothetical energy dispersion of Eq. 1|10|) and 
(b) that for doped SWNT. From this figure, it follows 
that the difference of ground state energies is a constant 
independent of iV$, because it corresponds to the sum of 
energy of all states for dispersion Eq. (|10|l within the re- 
gion of —it < X < ir and ir/2 < Y < ir. This shows that 
the persistent current in the doped armchair SWNTs has 
the same periodicity as the square-lattice cylinder after 
the replacement t — *■ V n , M — * 2n and N + 1 — > It', 
the period is thus given by &o/Z for an integer value of 
Z = Ix/n. Moreover, even when It In is equal to a ra- 
tional number Z /q for coprime integers Z and q, such 
as Ir/n = 3/5, the system shows an approximate frac- 
tional flux periodicity Qq/Z. It is because in Eqs. JSJ and 
(|f I \ Cz, C*2z, ■ • • are much larger than other Cj, owing 
to cos(27rZyl„) = I. 
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FIG. 2: Ground state energy as a function of the flux for (n, n) 
armchair SWNT. (a) n — 30 (diameter ~ 4 nm) with two 
different lengths: Ir/n = 3 (length ~ 22 nm) (solid curve), 
Ir/n = 9 (dashed curve), and (b) n — 80 (diameter ~ 10 nm) 
with two different lengths: It/ti = 0.6 (length ~ 11.6 nm) 
(solid curve), It /n = 1.2 (dashed curve). We offset the origin 
of the energy for comparison. 



with use of the result for the square lattice. Equation ||TJ 
for the square lattice is similar to Eq. © for the doped 
SWNT. Although there remains a factor of cos Y differ- 
ence, it does not affect the fractional periodicity; we can 
analytically show below that the fractional periodicity 
3>o/Z occurs when It / n = Z is an integer. To see this, 
it is useful to consider a hypothetical energy dispersion, 



-2I4| cos Y | (cos X + cos Y). 



(10) 



This energy eigenvalue corresponds to multiplying the 
factor \cosnn/{N + 1)| to the E nm (N&) for the square 
lattice. The corresponding persistent current can be ex- 
pressed by multiplying Cj in Eq. by the factor of 
|cosn?r/(iV + 1)|, 



C, 



N 

E 



N + 1 



sin(^)co><2,T,.l,,.,. (11) 



Thus the persistent current for Eq. l|lUf) still exhibits the 
fractional periodicity depending on the aspect ratio. To 
compare Eq. J^J with Eq. I|1(J|I . we consider the difference 



FIG. 3: Occupied states are shown in gray color in 
the Brillouin zone for the energy dispersion relations; 
(a) -2K|cosY|(cosX + cosK) in Eq. JTOJ and (b) 
-2Kcosr(cosX + cosY) in Eq. ©• 

Let us comment on the magnetic field B n for (n, n) 
armchair SWNT which corresponds to $o = 4 x 
10 -7 [gauss ■ cm 2 ]. For n = 100 (diameter is about 13 
nm), B n is about 30T, which is experimentally attain- 
able. For the fractional periodicity, the period becomes 
B n /Z ~ 30T/Z , which is easily achieved in experi- 
ments for larger Z. Although Z can be large in the 
long length limit of an armchair SWNT (It S> n), the 
periodic motion of the electron along the axis of long 
nanotubes may be affected by decoherence effect such 
as lattice deformations or defects. It is also valuable to 
comment on the magnetic moment (/i rbit) of the sys- 
tem^. The magnetic moment can be calculated directly 
from the persistent current as /i rbit = S , / pc (A r $), where 
S is the cross-sectional area of the tube (S = |C/ l | 2 /47r). 
When It/ti = Z (Z is an integer), we estimate that the 
maximum amplitude of the magnetic moment scales as 
Morbit/Ms = 0(l)n 2 , where ug = 6 x 10~ 2 [meV/T] is the 
Bohr magneton. 

The Fermi energy (E F = ±K - ±2.9eV) for doped 
SWNT for the fractional flux periodicity might be dif- 
ficult for chemical doping, whereas it might be possi- 
ble for electrochemical doping; a recent experiment of 
electro-chemical gating achieves Fermi energy shift of or- 
der ±leVi&. 
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Finally we point out that the fractional nature can 
be seen even for the smaller Fermi energy. In Fig. 01 
we plot Eo(N<s>) with a weak doping [Ep = — V^/10) 
for for (60, 30) (solid curve) and (70, 30) (dashed curve) 
chiral nanotube with It — 15 and It — 5, respectively. 
These curves clearly show some coherence effects. Thus a 
reproducible coherence pattern of the magnetic moment 
may be obtained as a function of magnetic field for doped 
SWNTs. 
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FIG. 4: Ground state energy of a weak doped (Ep = — K-/10) 
chiral nanotubes. We consider (60, 30) chiral nanotube with 
It = 15 (solid curve) and (70, 30) with It = 5 (dashed curve). 
The origin of the energy has been shifted to facilitate com- 
parison. 

It is not easy to obtain of a SWNT with the diame- 
ter of 4nm. Thus we generally need a multi-wall carbon 
nanotube in which the current flows only few outermost 



layers. In this case we need to consider an effect of a 
finite thickness of layers. If the periodicity comes from 
the some interlayer coupling, complex coherence effect 
would appear. In fact, there are several experimental re- 
ports on an approximately fractional flux periodicity in 
magneto-resistancei2*i2iiS. Bachtold et alii£ observed an 
oscillation with a period ~ <J>o/10. It not be explained by 
the Al'tshuler-Aronov-Spivak theory which predicts the 
period of $o/2 17,18 . At present moment, our results do 
not account for these oscillations, which will be a future 
work. 

In conclusion, we calculated numerically and analyt- 
ically the fractional periodicity of the ground state en- 
ergy and persistent currents in 2D cylinders composed 
of a square and doped armchair SWNT. A doped arm- 
chair SWNT also exhibits the fractional periodicity when 
they are doped to Ep = ±V n . The fraction ($ Q /Z) de- 
pends on the aspect ratio given by Z — It/™, where It 
is the number of unit cells along the cylindrical axis. An 
experimental investigation of the AB effect in a doped 
SWNT gives a key to understand this special coherence 
phenomenon. 
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